Astronomy & Astrophysics manuscript no. Outflow structure and velocity 'field of Orion Source!' new 
October 10, 2012 



©ESO 2012 



Outflow Structure and Velocity Field of Orion Source I: ALMA 
Imaging of SiO Isotopologue Maser and Thermal Emission 

F. Niederhofer, E.M.L. Humphreys, and C. Goddi 

European Southern Observatory, Karl-Schwarzschild-StraBe 2, D-85748 Garching bei Miinchen, Germany 
e-mail: fniederh@eso.org , ehumphre@eso.org, cgoddi@eso.org 



Preprint online version: October 10, 2012 



ABSTRACT 



Using Science Verification data fi-om the Atacama Large Millimeter/Submillimeter Array (ALMA), we have identified and imaged 
five rotational transitions (J=5-4 and J=6-5) of the three silicon monoxide isotopologues ^^SiO v=0, 1, 2 and ^^SiO v=0 and ^^Si^^O 
v=0 in the frequency range from 214 to 246 GHz towards the Orion BN/KL region. The emission of the ground-state ^^SiO, ^^SiO 
and ^^Si^^O shows an extended bipolar shape in the northeast- southwest direction at the position of Radio Source I, indicating that 
these isotopologues trace an outflow (~18 km s~\ P.A. ~50°, ~5000 AU in diameter) that is driven by this embedded high-mass 
young stellar object (YSO). Whereas on small scales (10-1000 AU) the outflow from Source I has a well-ordered spatial and velocity 
structure, as probed by Very Long Baseline Interferometry (VLBI) imaging of SiO masers, the large scales (500-5000 AU) probed 
by thermal SiO with ALMA reveal a complex structure and velocity field, most likely related to the effects of the environment of the 
BN/KL region on the outflow emanating from Source L 

The emission of the vibrationally-excited species peaks at the position of Source L This emission is compact and not resolved at an 
angular resolution of ~r/5 (~600 AU at a distance of 420 pc). 2-D Gaussian fitting to individual velocity channels locates emission 
peaks within radii of 100 AU, i.e. they trace the innermost part of the outflow. A narrow spectral profile and spatial distribution of 
the v=l J=5-4 line similar to the masing v=l J=l-0 transition, provide evidence for the most highly rotationally excited (frequency 
> 200 GHz) SiO maser emission associated with Source I known to date. The maser emission will enable studies of the Source I 
disk-outflow interface with future ALMA longest baselines. 

Key words. ISM: individual (Orion KL) - ISM: jets and outflows - stars: formation - masers 



1. Introduction 

Radio Source I in the Becklin-Neugebauer/Kleinmann-Low 
(BN/KL) region (d=418±6 pc; |Kim et "aLl|2008| ) provides the 
closest example of high-mass star formation. It is a deeply em- 
bedded young stellar object (YSO) and no infrared emission has 
been detected up to 22 fim ( Greenhill et al.||2Q04a ). However, 
Qkumura et al. ( 201 1^ detected a mid-infrared counterpart to 
Source I as a peak in the 7.8yum/12.4yum color temperature, in- 
dicating its protostellar nature. Measurements of the proper mo- 
tions of Source I and the BN object, the brightest mi d-infrared 
source with a bolometric luminosity of 1.3 x 10"^ Lq ( |de Buizer 
|et al.|2012] ) in this region, in combinati on with N-body simula- 
tions performed by|Goddi et al.|([2011a), suggest that Source I 
is a binary system with a total mass of -20 Mq that sustained a 
close (-50 AU) encounter with the BN object about 500 years 
ago. This dynamical interaction caused binary hardening (r < 
10 AU) while preserving the original circumbinary disk around 
Source I (| Moeckel"& Goddi„2012 ), and ejection of the high- 
energy, high- velocity, wide-angle outflow in the Orion BN/KL 
region observed in molecular hydrogen and carbon monoxide 
( |Bally et al.|2QTT] [Zapata et al.|2009| ). 

Source I powers a broad ensemble of SiO and H2O masers 
within radii between 10 and 1000 AU and is one of only three 
star-forming regions known to host SiO masers ( [Hasegawa et al. 
|1986 ). Vibrationally-excited v=l and v=2 SiO masers trace the 
innermost part of Source I inside 100 AU. Very Long Baseline 
Interferometry (VLBI) observations revealed an X- shaped pat- 
tern of the SiO masers that brackets Source I ( jGreenhill et al.| 



[T998l[Kim et al.|2008t [Matthews et al.|2010| ). In this structure, 
the north and the west arms are red- shifted, whereas the south 



and the east arms are blue-shifted. Additionally [Matthews et al. 
( 2010) found a bridge traced only by the v=2 masers that con- 
nects the south and the west arm and shows an intermediate ve- 
locity. The emission from the v= 1 and v=2 SiO is interpreted as 
arising from the limbs of a wide angle, funnel-like, slow (-18 
km s"^), bipolar outflow with a northea st- southwest axis that is 
launched from a rotatin g edge-on d i sk ([Greenhill et al. |2004bj 
Matthews et al. 2010). [Reid et al.[ ( [20Q7| ) observed a Al mm 
elongated continuum emission of -50 AU in radius that orig- 
inates from the center of the SiO maser emission, interpreted 
as an ionized disk elongated northwest-southeast. Monitoring of 
the Al mm continuum emission reveals a remarkably constant 
flux density and shape over a timescale of about 9 years, sup- 
porting the ionized disk model ( Goddi et al.|2011aj . 

Downstream the flow, at radii between 100-1000 AU, masers 
from H2O and vibrational groundstate SiO as well as thermal 
SiO emission trace a bip olar collimated, low density outflow 
originating from Source I ( Greenhill et al. 1998 2012 Plambeck 
et al.|2009| ). [Greenhill et al.[ ( [2012[ ) related this large-scale out- 



flow to a continuation of the wind emanating from the disk, 
making Source I an excellent case of a disk-wind collimated in 
a bipolar outflow. VLBI monitoring of SiO masers enabled for 
the first time a spatially resolved study of the region where the 
wind is launched. Although we have a relatively good knowl- 
edge about this inner region, the large scale structure of the 
outflow remains ambiguous. On scales > 1000 AU, the outflow 
shows a complex morphology and a velocity field that is incon- 
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sistent with the simple picture of a linear, collimated and rotat- 
ing outflow. Observations of the SiO v=0 J=2-l line performed 
with CARMA by [Plambeck et~aLl ([2009) show increasing ab- 
solute line-of-sight velocities with distance away from Source I 
and extended emission almost perpendicular to the outflow (cf. 
their Figure 3). [Plambeck et al.| ( |2009| ) interpreted this emission 
as due to a precession of the outflow on timescales of a few hun- 
dred years. This timescale is, however, inconsistent with a hard 
binary with an orbital separati on of a few AUs and an orbital pe- 
riod of a few years ( Goddi et al .|2011a Bally et al.|2011| ). Zapata 
et al.| (2012[) reported a subset of the ALMA data shown here, for 
the J=5-4 v=0 SiO line only, suggesting rotation of the outflow. 
The proposed rotation, however, is in the opposite sense to the 
rotation measured on smaller scales with SiO J=l-0 masers and 
which therefore challenges their interpretation. 

We report results from an ALMA Science Verification (SV) 
spectral survey at 214-246 GHz towards the Orion BN/KL re- 
gion. We identified and mapped 5 transitions of SiO isotopo- 
logues in the outflow of Source L Our main goal was to study in 
detail the large-scale structure and velocity field of the Source I 
outflow. These data also enabled us to find the first evidence for 
high-frequency (> 200 GHz) SiO maser emission at the position 
of Source L 

The structure of this paper is the following. In Section [2] the 
details of the observations as well as the imaging technique are 
described. The results are given in Section[3] In Section|4]we dis- 
cuss the structure and the velocity field of the outflow of Source 
I as well as the nature of the SiO emission detected with ALMA. 
Conclusions and future work prospects are drawn in Section [5] 

2. Observations and Imaging 

The ALMA SV spectral survey of the Orion KL Hot Core 
Region was conducted on 2012 January 20 and covered the 
frequency range 214-246 GHz which corresponds to the lower 
2/3 of the total frequency range in ALMA Band 6. 16 anten- 
nas (each 12 meters in diameter) were included in the compact 
Early Science (ES) Cycle configuration. Three observational 
runs were performed consecutively. The precipitable water vapor 
(PWV) values during the individual runs were ~2.3-2.8, ~1.8- 
2.6 and -1.8-2.5 mm respectively. Each observation consists of 
five diff'erent spectral setups including 4 spectral windows each. 
The resulting 214-246 GHz frequency range is covered by 20 
spectral windows of 1.875 GHz with an overlap of -0.2 GHz. 
Each spectral window consists of 3840 channels that are 488.281 
kHz wide. This corresponds to a Hanning smoothed resolution 
in velocity of -0.7 km s"^ The theoretical line sensitivity of the 
observations, estimated with the ALMA Sensitivity Calculator, 
is 0.01 Jy beam"^ in 0.7 km s"^ Table [l] summarizes the obser- 
vational parameters. 

The ALMA SV data were calibrated prior to public release. 
Callisto served as the absolute flux calibrator and the quasar 
J0607-085, with a flux of 1.4 Jy, was used as both bandpass and 
phase calibrator. For imaging we used the Common Astronomy 
Software Applications (CAS A). The first step of the SiO line 
imaging process was the continuum subtraction. We identified 
all channels with no line emission in each spectral window con- 
taining an SiO line and used the uvcontsub task which fits a 
model to the continuum and subtracts it afterwards. For the 
cleaning of the line datasets we used a pixel size of 0'.'3, an im- 
agesize of 30'' x 30'' and a Briggs weighting with a robustness 
parameter of 1. This provided a synthesized beam of -r.'7 x V!2 
at a P. A. of -1°. In each data cube, we used the strongest emis- 
sion feature to perform a single step, phase-only self-calibration. 



Tests with further steps of self-calibration did not improve the 
image quality significantly. The rms noise in the line-free chan- 
nels was about 0.01 Jy beam"^ while it was up to 0.3 Jy beam"^ 
in the channels with the strongest emission. 



3. Results 

We detected and imaged five lines in total from the three SiO 
isotopologues ^^SiO, ^^SiO and ^^Si^^O. All lines are J=5-4 ro- 
tational transitions except for the J=6-5 transition of the ^^Si^^O 
isotopologue. The ^^SiO and ^^Si^^O isotopologue transitions 
are in the vibrational ground-state whereas the ^^SiO is detected 
in the ground-state as well as in the excited states v=l and v=2. 
The parameters of these data are given in Table[2] Figure[T]shows 
the spectral profiles of the individual transitions, which appear 
in some cases to be blended with other molecular transitions. 
Complementary spatial information provided by the interfero- 
metric images in diff'erent velocity channels helped us distin- 
guish between diff'erent molecular lines and simplified the line 
identification. The "contaminants" are indicated in the spectrum 
with the name of the molecule at the respective line frequencies. 
The channel maps for each transition are given in the Figures 
[2] -[6] The channels are 2.5 km s"^ wide (velocity-averaged for 
imaging purposes) and cover the velocity range from -20 km s"^ 
to 40kms"^ 



-°SiO v=1 J = 5-4 





SiO v=0 J = 5-4 
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SiO v=0 J = 5-4 - 




'Si v = J = 6-5 




Velocity, Radio Defn LSR (km s ) 

Fig. 1. Spectra of the five J=5-4 (SiO) and J=6-5 (Si^^O) transi- 
tions from SiO isotopologue observed with ALMA. The dashed 
red line is at the systemic velocity of 5.5 km s"^ To guide the 
readers eyes, we overplotted the ^^SiO v=2 spectrum with three 
Gaussians of arbitrary width to highlight the three peaks that are 
due to the SiO emission. 
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Table 1. ALMA Observational Parameters 



ALMA Observational Parameters 



Number of Antennas 16 in compact ES Cycle configuration 

Frequency Range 214 -246 GHz (lower 2/3 of total frequency range in ALMA Band 6) 

Weather Conditions PWV (mm) -2.3-2.8 (Obs. Setup 1), -1.8-2.6 (Obs. Setup 2), -1.8-2.5 (Obs. Setup 3) 

Synthesized Beam - V.' 7x172 

Number of Spectral Windows 20 (3,840 channels each) 

Spectral Resolution -0.7 km s~^ 

Theoretical Sensitivity 0.01 Jy beam"^ in 0.7 km s~^ 



Table 2. Parameters of the Emission Lines of the isotopologues ^^SiO, ^^SiO and ^^Si^^O observed with ALMA 



Molecule 


Restfrequency 


Transition Synthesized Beam 


Peak flux 


Line Width « 


Brightness Temperature 




(GHz) 


J«-J/ OmC) X OmC); P.A.(°) 


(Jy) 


(kms-^) 


(K) 



29SiO V=0 


214.38574 


5-4 


1.76x1.18; 0.91 


33.46 


28 


>100 


SiO v=2 


214.08854 


5-4 


1.77x1.19; 0.46 


2.25 


29^ 




SiO v=l 


215.59595 


5-4 


1.75x1.18; 0.98 


44.97 


28 


>580 


SiO v=0 


217.10498 


5-4 


1.72x1.17; 1.29 


72.96 


55 


>50 


Si^^O v=0 


242.09493 


6-5 


1.63x1.08; 1.03 


5.42 







•^"^ Quoted is full width at zero maximum 
We estimated a lower limit of the velocity range, however due to line blending the line may in fact be broader 



3.1. The^^SiOv=1 Line 

The spectrum of the ^^SiO v=l J=5-4 line is clearly double 
peaked with the blue-shifted peak being stronger than the red- 
shifted one (Figure [T]). The velocity range of the emission is 
— 8 km s"^ to ~-f20 km s"^. The peaks themselves, located at 
km s"^ and 12.2 km s"^ respectively, are quite narrow with 
FWHM values of ~3.5 km s"^. The high velocity feature at the 
red side of the line emission is likely due to methyl-formate 
(CH3QCHQ). [Goddi et al . ( 2009a b) report observations of the 
J=l-0 line from SiO isotopologues towards Orion BN/KL at sev- 
eral epochs using the GBT and the VLA. For the masing ^^SiO 
v=l J=l-0 transition, they found a similar double peaked spec- 
tral profile at all epoch s with peaks at arou nd -5 km s"^ and 15 
km s"^ (cf. Figure 1 in [Goddi et al.|2009bt . The ALMA image 
shows that the v=l J=5-4 emission is peaked at the position of 
Source I and is not resolved with a beam of 1^7x1^2 (-700x500 
AU at 420 pc) at 215 GHz (Figure[2]). 

In order to characterize the structure and the kinematics of 
the ^^SiO v=l emission, we fitted a two-dimensional elliptical 
Gaussian model to the unresolved emission, in every velocity 
channel with emission greater that 4 times the rms noise in the 
peak line channel (i.e. 1.16 Jy beam"^), in at least six connected 
pixels. Since individual spots at similar velocities and diff'er- 
ent locations within the ALMA beam are blended together, the 
Gaussian fits give the intensity-weighted mean position of all 
the emission features in every velocity channel. Relative posi- 
tional errors are given by SO = 0.56b/S NR, where 6b is the 
FWHM of the synthesized beam. Errors range from 10 mas for 
the brightest channels to 50 mas in the channels with weakest 
emission. The fitted positions are shown in Figure [7] We see a 
clear velocity gradient in the northwest-southeast direction, with 
the blue-shifted emission being in the east and the red-shifted in 
the west. The same blended structure is seen in the v=l,2 ^^SiO 
and ^^SiO and ^^SiO v=0 J=l-0 transitions observed with the 
VLA by [Goddi et al.| ( [200 9b ). The velocity structure obtained 
from our Gaussian fits is consistent with the one found for the 
masing ^^SiO v=l and v=2 J=l-0 transitions by Matthews et 
al. (2010) despite the fact that we cannot resolve the detailed 
four-armed structure of the wind with our ALMA synthesized 
beam of VflxVf2. Nevertheless, the fitted positions are all lo- 
cated within ~100 AU from Source I. Hence, the vibrationally- 



excited v=l J=5-4 transition traces the innermost part of the out- 
flow associated with Source I, similarly to the v=l and v=2 J=l- 
transitions ( [Matthews et al.|2010| ). 

3.2. The^^SiOv=2Line 

The spectrum associated with the ^^SiO v=2 J=5-4 line is quite 
broad and shows several peaks (Figure |T]). However, only the 
three features at -6 km s"\ -hi 7 km s""and -h23 km s"^ are 
due to SiO emission. To guide the readers eyes, we have over- 
plotted the spectrum with three Gaussians of arbitrary width 
at the positions of each SiO line peak. The channel maps of 
the emission are given in Figure |3] The SiO v=2 emission is 
clearly unresolved, is peaked at the position of Source I, and 
probably also traces the innermost part of the outflow, similar 
to the J=l-0 v=2 line ( [Matthews et al.|2010| ). In some channels, 
we also see difl'use and quite strong emission that we attribute 
to ethyl-cyanide (CH3CH2CN), dimethyl-ether (CH3OCH3) and 
propanal (CH3CH2CHO). This emission is blended with the v=2 
SiO line and contributes to the remaining peaks in the spectrum. 
From the position of the peaks in the spectrum, we can set an 
interval from -6 km s"^ to -1-23 km s"^ as a lower limit to the 
velocity range of the ^^SiO v=2 emission. The spectrum o f the 
^^SiO v=2_J=l-0 transition obtained with the VLA by G oddi[ 
et al.[ ( [20Q9b| ), also shows a triple peaked profile, however with 



peaks at approximately -6 km s , -h6 km s and -h20 km s 
(cf. their Figure 1). 



3.3. The^^SiOv=0Line 

The ^^SiO v=0 J=5-4 line is the strongest and most extended 
line of all. In contrast to the other SiO lines, it does not have 
a pronounced double-peaked shape (Figure [T]). The peak of the 
emission is at about 6 km s"^ which is, with the given spectral 
resolution of 0.7 km s"^ consistent with the systemic velocity 
of Source I (-5.5 km s"^). The spectral profile shows an ad- 
ditional shoulder at 17 km s"^ The whole line emission is quite 
broad in velocity with a FWHM of about 22 km s"^ The channel 
maps for this line emission are presented in Figure 4] The ^^SiO 
v=0 J=5-4 line emission traces an elongated outflow structure 
in the northeast- southwest direction and peaks at the position 
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Fig. 2. Channel maps of the ^^SiO v=l J=5-4 line emission. The white contours are at 4, 8, 16, 32, 64 and 128 times the rms noise in 
the peak line channel (0.29 Jy beam"^). Every channel is 2.5 km s"^ wide and the central LSR velocity is indicated in the lower right 
corner of each box. The two narrow emission features in the spectral profile at around km s"^ and 12 km s"^ are clearly evident in 
the channel maps. The emission from the v=l ^^SiO line in the central channel maps is unresolved and centered at Source I (marked 
with a white star in each channel). The more diff'use emission at velocities > 25 km s"^ instead does not peak at the position of 
Source I and is likely associated with CH3OCHO. The ALMA synthesized beam is indicated as a white ellipse in the lower left 
panel. 



of Source I for velocities smaller than 20 km s"^. At velocities 
< 10 km s"^ the channel maps show an arc-like structure north 
of the position of Source I. In addition, there is also emission 
at velocities > 20 km s"^ that is not associated with Source I. 
The position of the bulk of this high- velocity emission consid- 
erably overlaps with the mid-infrared source IRc4. For analysis 
purposes, we created maps of the velocity-integrated emission 
(moment 0; Figure [8]) and the intensity- weighted velocity (mo- 
ment 1; Figure [9]). The moment map is obtained integrating 
the emission at each pixel in the velocity range from -20 to 40 
km s"\ while the moment 1 shows only pixels containing emis- 
sion that is above a threshold of 4 cr (= 1.08 Jy beam"^). 



3.4. The^^SiOv=0Line 



The spectrum of the ^^SiO v=0 J=5-4 transition also shows a 
double-peaked profile with two broad peaks at ~1 km s"^ and 
~14 km s"^ (Figure [1]). The emission shows a similar velocity 
extent to the ^^SiO v=0 J=5-4 transition. The blue end of the 
^^SiO spectrum appears to be contaminated by emission from 
the methyl-cyanide isotopologue ^^CHsCN originating from the 
hot core south of Source I. The channel maps (Figure [5| show 
an emission structure similar to the v=0 ^^SiO J=5-4 emission. 
As for the ^^SiO v=0 line, we created maps of the velocity- 
integrated emission and intensity- weighted velocity (Figures 10 
and The moment image shows the velocity-integrated 



emission from every pixel and the moment 1 map contains only 
pixels above 4 cr (= 0.68 Jy beam"^). 
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Fig. 3. Channel maps of the ^^SiO v=2 J=5-4 line emission. The contours in each channel are 3, 6, 9, 12, 15, 18, 21 and 24 times the 
rms noise in the peak line channel (0.065 Jy beam"^). Every channel is 2.5 km s"^ wide and the central LSR velocity is indicated in 
the lower right comer of each box. In the central channels, the SiO emission is unresolved and peaks at the position of Source I (pink 
star symbol), probably tracing the inner parts of the outflow associated with Source I (as the ^^SiO v=l transition). In addition, strong 
diflTuse emission is evident away from the systemic velocity blending the SiO line, that we attribute to CH3CH2CN, CH3CH2CHO 
and CH3OCH3. The ALMA synthesized beam is indicated as a white ellipse in the lower left panel. 



3.5. The^^Si^^Ov=0Line 



4. Discussion 



The ^^Si^^O spectral profile shows a double-peaked line shape 
similar to the other isotopologues. The blue-shifted peak, how- 
ever, appears to be highly blended by CH3CH2CN emission orig- 
inating from the Hot Core (see for example the channel maps in 
Figure [6]). We cannot determine the exact velocity range of the 
SiO emission, because of the strong blending in the blue part of 
the emission. Nevertheless, from the channel maps one can es- 
timate that the SiO emission ranges from about 2 km s"^ up to 
17 km s"^ In this velocity range, the ^^Si^^O emission shows a 
structure that is similar to, but weaker than the ^^SiO line emis- 
sion. 



4.1. The Structure and Velocity Field of the Outflow 

The outflow from Source I has a well-ordered spatial and ve- 
locity structure within 1000 AU, as evidenced by high- angular 
resolution imaging of the SiO J=l-0 rotational lines in the 



vibrationally-excited states (Matthews et al. 20101 and the 
ground-state (Greenhill et al.||2012|). In particular, the flow 



appears as a wide-angle rotating wind within 100 AU from 
Source I, and then collimates into a bipolar, linear, rotating flow 
at radii of 100-1000 AU. The new ALMA images of the J=5- 
4 rotational lines from diff'erent isotopologues reported here, 
probe larger scales of the flow, up to -5000 AU. Interestingly, 
these large scales show a rather complex outflow morphology 
and velocity field when compared with the innermost flow re- 
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Fig. 4. Channel maps of the ^^SiO v=0 J=5-4 line emission. The contours in each channel are 3, 6, 9, 12, 15, 18, 21, 24, 27, 30 and 
33 times the rms noise in the peak line channel (0.27 Jy beam"^). Every channel is 2.5 km s"^ wide and the central LSR velocity 
is indicated in the lower right corner of each box. The bulk of the emission comes from an extended structure elongated northeast- 
southwest and centered at Source I (white star in each channel). The ALMA synthesized beam is indicated as a white ellipse in the 
lower left panel. 



gions (Figures [8] |9][T0|[TT]). We interpret this peculiar behaviour 
on different scales as a result of the interaction of the outflow 
from Source I with the complex environment of BN/KL. In the 
following subsections, we will focus on investigating the rela- 
tion of the ground- state ^^SiO and ^^SiO emission probing the 
large- scales of the outflow with the infrared sources located in 
the vicinity of Source I, the explosive event occurred about 500 
years in BN/KL, and dense gas in the Orion hot core. 

4.1 .1 . The v=0 ^^SiO Emission 

Figures [8] and [9] show the moment and moment 1 maps of the 
^^SiO J=5-4 transition observed with ALMA, overlaid with the 
contours of the total intensity of the ^^SiO v=0 J=2-l emission 
as observed by Plambeck et aL] ( |2009| ) using CARMA. The 5- 
4 and 2-1 transitions of ^^SiO are broadly consistent both in 
morphology and position angle, and clearly trace a bipolar out- 



flow driven by Source I into its surrounding molecular cloud 
along a northeast- southwest axis. Although the outflow has a 
pronounced bow-tie morphology, it is not symmetric. The south- 
ern lobe appears broader and more extended than the northern 
one, whereas, perpendicularly to the outflow axis, the emission 
in both lobes is more extended to the north-west (Figure [5]). 

The velocity field in the moment 1 map displays more inter- 
esting details on the large-scale structure of the outflow (Figure 
|9]). In particular, three new elements are observed. 

i) The southern lobe displays two wings of red-shifted emis- 
sion, the first one to the east is seen in the velocity range 
from ~7.5 km s"^ to ~ 17.5 km s"^ and the second one to the 
west traces the highest velocity red-shifted gas, from ~20 
km s"^ to ~35 km s"^ ("high-velocity wing"; Figure|9]). 

ii) The end of the north lobe of the outflow appears to be "bent" 
to the east resulting in an arm-like structure of faint highly 
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Fig. 5. Channel maps of the ^^SiO v==0 J=5-4 line emission. The contours in each channel are 4, 8, 16, 32 and 64 times the rms noise 
in the peak line channel (0.1 Jy beam"^). Every channel is 2.5 km s"^ wide and the central LSR velocity is indicated in the lower 
right corner of each box. The emission peaks at Source I (white star) and shows a bipolar structure elongated northeast- southwest. 
It is double-peaked in velocity space with maxima at 1 km s"^ and 13 km s"^ The ALMA synthesized beam is indicated as a white 
ellipse in the lower left panel. 



red-shifted (>25 km s"^) emission, having a northwest- 
southeast orientation ("red arm"; Figure [9]). 
iii) An arc-like feature of blue-shifted emission ( — 10 km s"^ 
to ~5 km s"^), centered on Source I, extends along the 
northwest-facing edge of the outflow ("blue arc"; Figure|9]). 

Since these features are evident in other t hermal transitions 
of SiO, like thev=0J=2-l ( |Plambeck et al. ^2009 ) and the v=0 
J=8-7 (Zapata et al.|2012 ), we exclude contamination from other 
molecular species, i.e. they must be associated with SiO emis- 
sion. 

Consistently with our findings, [Zapata et al.| ( |2012| ) noticed 
that the J=5-4 v=0 SiO emission has a butterfly-like shape with 
the blue-shifted emission in the northwest-facing side of the out- 
flow (corresponding to our blue arc) and the red-shifted emission 
in the southeast-facing side (corresponding to our moderately 
high- velocity wing); they also noticed a high- velocity "tail" ly- 
ing in between the two wings (see their Figure 3), which cor- 
responds to our "high- velocity wing" (Figure [9]). They discuss 



some possibilities to explain this peculiar velocity structure, in 
particular the simultaneous presence of two bipolar outflows 
from a binary or a single rotating outflow. Two outflows hav- 
ing north- south and east- west axes, respectively, could explain 
the blue and red wings but can be ruled out by the presence of 
the highest red- shifted tail or wing in the southern lobe. The 
second option of a single outflow assumes a sense of rotation 
with redshifted emission in the southeast-facing side and blue- 
shifted emission in the northwest-facing side of the outflow. On 
scales smaller than 100 AU, however, Ma tthews et a l. ( 2010) ob- 
served rotation in the opposite direction, requiring a change of 
the rotation sense from small to large scales, which is unlikely. 
Precession in a binary over a timescales of several hundreds of 
years could in principle explain the difl'erence in the red- shifted 
and blue- shifted velocities at difl'erent scales ( [Plambeck et al.| 
2009 ), but it is inconsistent with the strong evidence that Source 
I is a hard binary with an orbital separation of (at most) a few AU 
and an orbital period of a few years ( [Goddi et al.|201 la[|Bally et 
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Right Ascension (J2000) 

Fig. 6. Channel maps of the ^^Si^^O v=0 J=6-5 line emission. The contours are at 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44 and 48 
times the rms noise in the peak line channel (0.07 Jy beam"^). Every channel is 2.5 km s"^ wide and the central LSR velocity is 
indicated in the lower right corner of each box. The emission from the ^^Si^^O isotopologue is mostly compact and peaks on Source 
I (pink star), but displays also a low-brightness elongated structure in the northeast- southwest direction, indicating that the ^^Si^^O 
isotopologue traces the same outflow structures as the ^^SiO v=0 and ^^SiO v=0 lines. The diflTuse and extended emission at larger 
velocities mostly originates from the hot core, south of Source I, and may be attributed to CH3CH2CN. The ALMA synthesized 
beam is indicated as a white ellipse in the lower left panel. 



al. 2011| ). Besides, while we see a velocity oflfset from red to blue 
in the northern lobe in the southeast-northwest direction, in the 
southern lobe we observe the opposite trend, with the more red- 
shifted emission in the north-west ("high- velocity wing"; Figure 
|9]). The outflow rotation scenario faces a similar challenge in ex- 
plaining the tail with the highest red-shifted emission lying in 
between the mod erately red-shifted and the blue-shifted wings 
( [Zapata et al. (2012; their Figure 3). We conclude that on scales 
of thousands of AU there is not evidence for a clear velocity gra- 
dient perpendicular to the outflow axis that can be interpreted as 
rotation. 

We propose an alternative explanation, that the peculiar fea- 
tures observed in the SiO velocity field on large scales are most 
likely related to the outflow interaction with its environment. For 



this purpose, we overlaid the SiO velocity field map with mid- 
infrared sources in the vicinity of Source I (Figures [9] and \T3) 
These overlays evidence three interesting elements. 



i) The contours of the southern lobe of the SiO emission 
clearly follow the 12 yum emission from IRc4 (Figure p3]) and 
the high-velocity red wing significantly overlaps with IRc4 
(Figure [9]). Using SOFIA, de Buizer et al.| ( |2012| ) propose 
from their IR observations of the BN/KL region that IRc4 is 
self-luminous. This high red- shifted wing of SiO emission 
may not be associated with the outflow of Source I but with 
IRc4, that harbors a YSO that could heat the dust in its vicin- 
ity and enrich SiO via sublimation. Alternatively, the inter- 
action of the Source I outflow with IRc4 could create shocks 
that sputter dust grains enriching gas phase SiO. 
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20 kms 
> 20 kms"^ 



S^'as'" 14^520 



Fig. 7. Fitted positions of the ^^SiO v=l emission. The crosses 
are color coded by velocity. We see that there is a clear veloc- 
ity gradient along the northwest- southeast direction, with blue- 
shifted emission in the eastern part and red-shifted emission in 
the weste rn part. This is in good agreement with the result ob- 
tained by [Matthews et al.| ( [5UTQ| from VLB A data for the v=l 
J=l-0 transition. The size of the image is about 160 AU x 160 
AU. 




Fig. 8. Total intensity (moment 0) map of the v=0 J=5-4 ^^SiO 
line. Qverplo tted are contours of the v=0 J=2-l ^^SiO total in- 
tensity from [Plambeck et al.| ( |20Q9| ) obtained with CARMA. 
Contours are at 0.1, 0.2, 0.4, 0.7, 1.1, 1.5, 1.9, 2.6 and 3.7 
Jy beam"^ km s"^ . The white star marks the position of Source I. 
The axis of the outflow (P.A. -50°) is indicated by a white line. 



ii) The northern part of the bended high- velocity "red arm" aris- 
ing from the north lobe to the east of Source I overlaps well 
with the positions of the mid-infrared sources IRcll and 
IRcl2, which may play a role in exciting this high- velocity 
SiO emission. IRcll and IRcl2 alone, however, cannot ex- 
plain the whole "red arm". Since the latter has a filamentary 
structure with a northwest-southeast orientation, reminiscent 
of the "fingers" of shocked H2 bullets in the fast OMC-1 
outflow, we speculate that the "red arm" (in particular its 
southernmost part) may be tracing part of the high- velocity 




5 35 15 



Fig. 9. Intensity-weighted velocity (moment 1) map of the v=0 
J=5-4 ^^SiO line emission w ith overlaid contours fr om the v=0 
J=2-l 2^SiO total intensity ( [Plambeck et al.||2009| ). Only pix- 
els with fluxes >1 Jy beam^ or 4 times the rms noise of the 
channel with the brightest emission, were used. The white el- 
lipses indicate the position s ([Shuping et al.||2004 ) and approxi- 
mate sizes (estimated from ] Shuping et al.|2Q04 ) of the infrared 
sources IRc4, IRc7, IRcll and IRcl2. Source I is marked with 
a white cross. The newly discovered infrared source SOFl pe| 
,Buizer et al.,2012) is indicated as a circle in the upper left cor- 
ner. The two white crosses in the blue arc labeled 'Goddi' and 
'Zapata' are t he points of interacti on between Source I and BN 
as derived by [Goddi etaT] ( |20 1 1 at and Zapata et aL]p009) , re- 
spectively. 



northwest- southeast outflow that originated from a close en- 
counter o f Source I and the BN object approx imately 500 
years ago ( [Goddi et al.|2011a[[Bairy et al.|2011| ). 
iii) The center of the "blue arc" to the north-west of Source I 
falls in between two independent measurements of the cen- 
ter of dynamical interaction between Source I and the BN 
object ( [Goddi et al.][2011at [Zapata et al.][2Q0 9i Figure [9]). 
We speculate that this blue-shifted emission, similarly to the 
red- shifted emission from the "red arm", could arise from 
remnant material that was stripped away from the interact- 
ing objects with moderately velocity. 



4.1 .2. The v=0 ^^SiO Emission 



Figures 10 and 11 show the ALMA moment and moment 1 
maps of the ^^SiO J=5-4 transition, overlaid with the contours 
of the total int ensity of the CARMA ^^SiO v=0 J=2-l emission 
( Plambeck et al..2009j . The struc ture of ^^SiO J^5-4 lin e fol- 



Plambeck et al. 



(2009) quite 



lows the^^iO J=2-l contours by 

well (Figure 10). The ^^SiO J=5-4 emission appears more com- 
pact than the^SiO J=5-4 emission but has an overall similar 
morphology, with the southern lobe of the outflow being more 
extended than the northern lobe (e.g. compare Figures [8] and 10). 

The velocity field shows a much clearer blue-red asym- 
metry than the ^^SiO perpendicular to the outflow axis (along 
northwest- southeast). This apparent velocity gradient could be 
in principle interpreted as rotation. We argue, however, that this 
asymmetry, as well as the one displayed in the ^^SiO emission. 
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14 5 

RA (J2000) 



Fig. 10. Moment map of the v=0 J=5-4 ^^SiO transition over- 
plotted with contours of the v=0 J=2-l ^^SiO total intensity. The 
contour levels are the same as in Figure [S] The star symbol indi- 
cates the position of Source 1. 
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RA (J2000) 



Fig. 11. Moment 1 map of the v=0 J=5-4 ^^SiO line with over- 
laid contours from the v=0 J=2-l ^^SiO total intensity. Only 
pixels with fluxes larger than 4 times the rms noise of 0.17 
Jy beam"^ were used to produce this map. The crosses mark the 
positions of IRcll, Source I, IRc7 and IRc4. 



cannot be attributed to rotation but most probably to the interac- 
tion of the outflow with its environment. In particular, the blue- 
shifted emission along the northwest-facing edge of the outflow 
is observed also in the ^^SiO emission, for which we favor an 
association with the explosive event rather than Source I itself. 
On the contrary, the red-shifted emission (with velocities > 20 
km s"^) is mainly composed of two red features t hat are not 



Plambeck 



seen in the ^^SiO J=5-4 line or in the J=2-l line by 
et al.| ([2009|: 1) an elongated structure in the southeast-facing 



edge of the outflow; 2) a red clump near the position of IRc7. As 
this emission is not seen in the other SiO transitions, it may be 
associated with another contaminant molecule rather than SiO. 
Figure 12 shows the moment 1 map of the ^^SiO J=5-4 transi- 




ts 

RA (J2000) 



Fig. 12. Moment 1 map of the v=0 J=5-4 ^^SiO line overlaid 
with contours of the ammonia column density map from EVLA 
observations ( |Goddi et al.||2011b| ). The contours are from 35% 
to 90% of the peak column density of 2.0x10^^ cm"^, in steps of 
5%. The star symbol indicates the position of Source I and the 
cross marks the position of IRc7. 




Fig. 13. Orion BN/KL 12.5 jam mosaic ( Greenhill et al.|2004a ) 

overlaid with v=0 J=5-4 ^^SiO integrated intensity contours. The 
contours are from 10% to 90% of the peak emission of 170 
Jy beam"^ in steps of 5%. The crosses indicate the positions of 
Source I, Source n, IRc7 and the BN object (from left to right). 



tion overlaid with contours of the ammonia column density map 
( |Goddi et al.||2011b ) tracing dense gas in the hot core region. 
The contours appear to overlap significantly with the two red 
features. 

Since the high red- shifted ^^SiO emission seems to follow 
the dense gas traced by the ammonia column density (Figure 
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12 ), we suggest that it may trace the interface of the Source I 
outflow with the hot core and IRc7, resp ectively. 

29SiO 



Tercero et al. 



(2011) showed that the 



Consistently, 

v=0 J=5-4 line is contaminated by the following molecules: the 
hydrogen-sulfide isotopologue ^^SH2, the methy l-cyanide iso- 



topo logue ^^CHsCN, and sulfur-monoxide (SO). [Beuther et al. 
( [2QQ5| ) imaged the (196,i9-186,i8) ^^CHsCN line at 339,137 GHz, 
which clearly peaks both at the hot core and IRc7 positions. If 
the ^^CHsCN emission is responsible for the strong transition 
at 214.37437 GHz, we estimate an LSR velocity of > 6 km s"^ 
which is reasonable for emission coming from the hot core (v~9 
km s"^). We therefore propose that the apparent high- velocity 
features in the ^^SiO maps can be attributed to ^^CHsCN and be 
associated with IRc7 and the interface between the hot core and 
the outflow from Source I. 

We conclude that the blue-red asymmetry in the SiO moment 
1 map is only apparent. Therefore, we exclude the presence of 
a velocity gradient perpendicular to the outflow axis as well as 
evidence of outflow rotation on scales of thousands of AU from 
Source I. 



4.2. The Nature of the Emission 

In order to investigate the nature of the emission (thermal vs. 
maser), we estimated the brightness temperature of the ^^SiO 
v=0, v=l and the ^^SiO v=0 transitions (Table|2]). is given by: 
Th = (S A^)/(2kBdQ), where S is the flux, A the wavelength and 
dQ. the solid angle of the source. The ^^SiO v=l emission is not 
resolved, therefore we can only give a lower limit to the bright- 
ness temperature. For a beamsize of r.'7xlV2 and a peak flux 
of 45 Jy, we derive a brightness temperature -580 K. Based 
on the lower limit of the brightness temperature alone, the nature 
of the v=l emission remains ambiguous. However, the two spec- 
tral peaks of the ^^SiO v=l transition are narrow (line-widths 3.5 
km s"^ Figure[T]), consistent with maser emission. Moreover, the 
positions fitted to the unresolved emission at diff'erent velocities 
show a distribution similar to the J=l-0 v=l masing transition 
( [Matthews et al.|[2QlQ| ) (see Section 3.1 and Figure 17]). These 
elements provide first evidence that the v=l J=5-4 ^^iO tran- 
sition is a maser. If confirmed, this would be the first detection 
of high-frequency SiO maser emission (> 200 GHz) associated 
with Source I. In addition, the masing nature of the v=l J=5-4 
transition makes it a valuable tracer of the inner (< 100 AU) part 
of the Source I outflow for follow up observations using ALMA 
long baselines (15 km). 

The v=0 lines from ^^SiO and ^^SiO are resolved, so we cal- 
culated Tij by integrating the flux over the angular size of the 
source. We obtained brightness temperatures of ~50 K and -100 
K for the ^^SiO and the ^^SiO isotopologue, respectively. The 
low brightness temperature of the v=0 ^^SiO and ^^SiO transi- 
tions indicates that the emission is most likely thermal. But the 
fact that the emission from the ^^SiO isotopologue is brighter 
than the one arising from the ^^SiO molecule, suggests that the 
^^SiO transition could be partly masing. 



5. Conclusions and Future Work 

We have reported ALMA observations of five rotational tran- 
sitions of three SiO isotopologues towards the Orion BN/KL 
region. The vibrational ground- state emission from ^^SiO and 
^^SiO shows an extended bipolar structure along a northeast- 
southwest axis at the position of Source I and traces the outer 
parts of its outflow up to scales of a few thousand AU. These 



large- scales reveal complex morphological and dynamical struc- 
tures that may be due to the outflow interaction with its envi- 
ronment, in particular the hot core and the mid-infrared sources 
in the vicinity. Emission from the vibrationally-excited v= 1 and 
v=2 J=5-4 ^^SiO is not resolved with the current resolution 
(~ir5; -600 AU at 420 pc). 2-D Gaussian fitting to the un- 
resolved emission in each velocity channel indicates it traces 
the inner regions (<100 AU) of the outflow associated with 
Source I. Based on the narrow line- widths of spectral features 
(<3.5 km s"^) and a distribution similar to the masing ^^SiO 
J=l-0 v=l line, we found first evidence for the ^^SiO v=l J=5- 
4 transition being masing. Follow-up observations with the full 
ALMA array will enable resolution of the region of the ^^SiO 
v=l J=5-4 emission and provide data complementary to those 
taken at lower frequencies, possibly unveiling new dynamical 
and physical properties of molecular gas in a disk/outflow sys- 
tem associated with the closest known high-mass YSO. 
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